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How does a chemical reaction proceed at ultralow temperatiires? Can sim- 
ple quantiim mechanical rules such as quantiim statistics, single scattering 
partial waves, and quantum threshold laws provide a clear understanding for 
the molecular reactivity under a vanishing collision energy? Starting with an 
optically trapped near quantum degenerate gas of polar ^°K^^Rb molecules 
prepared in their absolute ground state, we report experimental evidence for 
exothermic atom-exchange chemical reactions. When these fermionic molecules 
are prepared in a single quantum state at a temperature of a few hundreds 
of nanoKelvins, we observe p-wave-dominated quantum threshold collisions 
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arising from tunneling through an angular momentum barrier followed by a 
near-unity probability short-range chemical reaction. When these molecules 
are prepared in two different internal states or when molecules and atoms are 
brought together, the reaction rates are enhanced by a factor of 10 to 100 due 
to s-wave scattering, which does not have a centrifugal barrier. The measured 
rates agree with predicted universal loss rates related to the two-body van der 
Waals length. 

Scientific interests in precisely understanding the fundamental aspects of chemical reactions 
and controlling their dynamic processes have stimulated pioneering work on molecular beams to 
study state-to- state reactions using molecular alignment, velocity selections, and angle-resolved 
measurement (1-4). However, substantial motional energies remained in earlier work and ther- 
mal statistical averages were a necessary ingredient. By preparing a molecular ensemble in the 
quantum regime, we expect to develop fundamental insights into how chemical reaction pro- 
cesses may be precisely guided by quantum mechanics. Reaction dynamics at vanishingly low 
energies remain a fascinating and yet unexplored scientific realm (5). Here both the internal and 
external molecular degrees of freedom are prepared quantum mechanically, and each step of a 
complex reaction may be analyzed based on single quantum states and single reaction channels. 

That chemical reactions could occur at ultralow temperatures seems at first glance counter- 
intuitive. However, ultracold collisions, where particles scatter only in the lowest angular mo- 
mentum partial wave, are governed by quantum statistics and quantum threshold behaviors 
described by the Bethe-Wigner laws (6-8). In this regime, particles are represented by their de 
Broglie wavelength, which increases with reduced temperatures. This wave nature of particles 
replaces our intuitive and classical picture of colUsions. The wave manifestation of tunneling 
through reaction or angular momentum barriers may play a dominant role in dynamics and scat- 
tering resonances can have dramatic effects on reactions (9). In addition, any barrierless chem- 
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ical reactions will always take place when two reactants are sufficiently close together {10). In 
this case, chemical reaction rates will be determined to a large extent by coUisional properties at 
large intermolecular separations, and thus by how the two partners approach each other. Once 
their separation reaches a characteristic length scale (~10 ao, with ao = 0.53 • 10"^° m), a 
chemical reaction happens with a near unity probability. Therefore, chemical reactions can be 
surprisingly efficient even at ultracold temperatures. Indeed, this model for barrierless reactions 
predicts loss rates that are universal in the sense that they do not depend on the details of the 
short-range interactions, but instead can be estimated using only knowledge of the long-range 
interactions (77). 

Like the case of collisions of ultracold atoms, the study of ultracold chemical reactions will 
play a fundamental role in advancing the field of molecular quantum gases. For example, under- 
standing and manipulating collisions of atoms at ultralow temperatures (< 1/iK) have been cru- 
cial for the realization of quantum degenerate gases (12-14), Fermi superfluids that provide op- 
portunities to explore the underlying connection between superconductivity and Bose-Einstein 
condensation (75), neutral atom-based systems for quantum information science (16-18), and 
strongly correlated quantum gases (79, 20). Ultracold molecules undergo a more diverse set 
of coUisional processes, with distinct inelastic collision mechanisms arising from chemical re- 
actions, in addition to the traditional state -changing collisions seen with ultracold atoms and 
highly vibrationally excited molecules (27). Furthermore, polar molecules possess anisotropic 
and long-range dipolar interactions that can be precisely controlled with external electric fields, 
with rich prospects of coUisional resonances (22, 23). Experimental investigation of molecu- 
lar collisions is essential for future applications including studying anisotropic and collective 
behavior in quantum gases (24, 25), modeling new quantum phases and exotic many-body 
physics (26), implementing schemes for quantum information (27), and developing tools for 
precision physical and chemical measurements. 
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We focus here on the study of ultracold collisions, including chemical reactions, of °K^'Rb 
molecules, which we prepare in their lowest electronic, vibrational, rotational, and hyperfine 
energy state at a high phase-space density {28, 29). We find clear evidence of the essential role 
that quantum statistics and quantum threshold laws play in determining the rates of inelastic 
collisions. Our experimental observations confirm a universal loss mechanism. 

A prerequisite for exploring ultracold chemical reactions is a gaseous molecular sample that 
is sufficiently dense, ultracold, and suitable for precise control of specific quantum states (5). 
The starting point for this work is an ultracold trapped gas of fermionic ''^K^^Rb molecules pre- 
pared in a single hyperfine level of the ro-vibronic ground state (A^ = 0, f = OofX^S+) {28, 
29). The optical trap depth is ~ ks-lO yuK, where is Boltzmann's constant. The molecules 
are produced using a single step of two-photon Raman transfer from extremely weakly bound 
molecules at a magnetic field of ~545.9 G. The coherent transfer is efficient and does not heat 
the gas, resulting in a gas of ro-vibronic ground-state molecules with an average number den- 
sity of 10^^ to 10^^ cm^"^ and a translational temperature of a few hundreds of nanoKelvin. At 
this ultralow temperature, even the tiny molecular hyperfine-state energy splittings are much 
larger than the translational energy. Manipulation of the hyperfine states hence becomes ex- 
tremely important and relevant in exploring possible collision channels. In addition, complete 
control over the internal quantum state of the molecules permits direct observation of the role 
of quantum statistics in determining the molecular interactions. 

A precise measurement of the hyperfine structure and the manipulation of individual hyper- 
fine state populations in the X^S+ ground state of '^*^K®''Rb were reported in (29). The X^S+ 
state has zero total electronic angular momentum, so that the hyperfine structure is basically 
the Zeeman effect of the nuclear spins /^^=4 and 1^^=312 {11, 30) at the applied magnetic field. 
The hyperfine structure is depicted in Fig. [T| where a total of 36 states are labeled by their 
projections of the individual nuclear spins, mf^ and mf. For the current study, we produce 
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molecules either in a single spin state (\mf ,mf^)) or in a mixture of two spin states. The hy- 
perfine states used are an excited state | — 4, 1/2) and the lowest energy spin state | — 4, 3/2); 
these are marked by the two ellipses in Fig. [Tj The | — 4, 1/2) state is populated directly by the 
two-photon Raman transfer starting from the weakly bound molecules {28, 29). The spin state 
of these molecules can be further manipulated coherently. For example, the entire | — 4, 1/2) 
population can be transferred into the lowest hyperfine state, | — 4, 3/2), using two successive 
TT— pulses through a rotationally excited = 1 intermediate level. The = 1 level has strong 
nuclear electric quadrupole interactions that couple rotations with nuclear spins, enabhng nu- 
clear spin flips (29, 30). We can probe molecules in any particular hyperfine state by reversing 
the entire transfer process and putting the population back into the initial weakly bound state. 
We then use high signal-to-noise-ratio absorption imaging to measure the molecular gas number 
and temperature. 

The ability to control molecular internal states including the electronic, vibrational, rota- 
tional, and nuclear spin degrees of freedom, and in particular the possibility of preparing them 
in the lowest energy state, allows us to probe inelastic collisions in a way that limits unwanted 
loss mechanisms and reduces ambiguities in the identification of possible chemical reaction 
channels. This is crucial because it is difficult in the ultracold gas experiment to find probes 
for the direct observation of reaction products. In Table [T] we consider the binding energies for 
various types of molecules made from different combinations of and ^^Rb atoms. These 
energy estimates allow us to assess whether a specific two-body reaction process is endothermic 
or exothermic. We note that all the results reported here were obtained in the absence of any 
external electric field, and hence the effective molecular dipole moment in the lab frame is zero. 

To probe the quantum nature of ultracold molecular collisions and chemical reactions, we 
begin by preparing the KRb molecules in a single nuclear spin state of the ro-vibronic ground 
state. All unpaired atoms that remain after the initial stage of the molecular creation process are 
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Table 1: Summary of the relevant molecular energetics involved in possible chemical reac- 
tions. The binding energies are given with respect to the threshold energy for free atoms in the 
absence of a magnetic field. The ^''Rb2 and '^°K2 binding energies include isotope shifts from 
the data in the respective references. The trimer binding energies are unknown. Calculations of 
trimer binding energies are needed and will be important for future experiments on any bi-alkali 
species. 



Molecules 


V = binding energy (Dq) 


reference 


''Rb2 


3965.8(4) cm-i 


(31) 


40K87Rb 


4180.417 cm-i 


(28) 


^°K2 


4405.389(4) cm-^ 


(32) 


K2Rb 


unknown 




KRb2 


unknown 





selectively removed from the optical trap via resonant light scattering (33). From an argument 
based on the energetics summarized in Table[T} the molecule-molecule collisions have a possible 
exothermic chemical reaction, namely KRb -i- KRb— i>K2 -i- Rb2, which releases ~10 cm^ of 
kinetic energy. The reactions KRb -i- KRb^K2Rb -i- Rb and KRb -i- KRb^KRb2 -i- K could also 
be exothermic. All of these reactions require breaking and making molecular bonds. If the KRb 
molecules are prepared in an excited hyperfine state, spin relaxation to a lower hyperfine state 
provides an additional inelastic scattering mechanism. 

The quantum statistics of the molecules plays an essential role in collisions at a temperature 
of a few hundred nanoKelvin, where collisions with large-impact parameters and correspond- 
ingly large centrifugal barriers are frozen out and the collisions are typically dominated by a 
single partial wave with orbital angular momentum quantum number L=0 (s-wave) or L=l (p- 
wave). Our KRb molecules are fermions and therefore the total wave function describing a 
KRb-i-KRb collision is antisymmetric with respect to molecular exchange. For spin-polarized 
molecules all prepared in exactly the same internal quantum state, the p-wave is the lowest 
energy symmetry-allowed collision channel. The height of the centrifugal barrier for the L=\ 
KRb-KRb collisions is A;^ ■ 24 fiK (34, 35). This barrier height is over an order of magnitude 
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larger than kBT, where T is the translational temperature of the molecular gas. Thus, collisions 
of spin-polarized molecules are expected to proceed predominately via tunneling through the 
p-wave barrier. Note that T is greater than 1.4 times the Fermi temperature. If two molecules 
make it through the barrier to short range, chemical reactions or hyperfine state-changing colli- 
sions can take place, leading to a loss of the entrance channel population. We note that even a 
single nuclear spin flip corresponds to a released quantity of energy that is above the trap depth, 
and would contribute to loss of trapped molecules. 

The quantum nature of the collisions can be seen in the temperature dependence of loss rates. 
The Bethe-Wigner threshold law predicts that the p-wave inelastic/reactive collision rate should 
be linear in temperature (oc T). To look for this behavior, we first prepared spin-polarized 
molecules in the single hyperfine state | — 4, 1/2) for various T ranging from 200 to 900 nK (38). 
The temperature is measured from the expansion energy of the molecular gas after releasing it 
from the optical trap. For each initial temperature, we observed the time-dependent molecular 
loss (39) and extracted a two-body loss rate P (which is 2 times the coUisional event rate) by 
fitting the measured decay of the molecular gas density n vs. time t (Fig. [2]^A)) to 

— = —on — an. (1) 

dt 

Here, the first term accounts for number loss and the measured P can be compared to theoretical 
predictions. The second term accounts for density change due to heating of the trapped gas 
during the measurement. Within a single measurement, we observe an increase in temperature 
that is at most 30%. In the analysis for each data set, we fit the measured temperature to a 
linear heating rate and obtain a constant slope c. In Eq. 1, we then use a = f^^^, where 
T is the initial temperature. At our lowest temperature of 250 nK, the heating was 7(1) nK/s 
and (3 = 3.3(7) • lO^^^cm^/s (Fig.jlj^A)). The measured dependence of /3 vs. T is summarized 
in Fig. ^B) (closed circles). Here, we fit the data to a power law /5(T) oc and find that 
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L = 1.1(2), which agrees with the predicted p- wave threshold law. This result demonstrates 
that indistinguishable ''°K^''Rb molecules at ultralow temperatures collide via tunneling through 
a p-wave barrier followed by an inelastic collision in the short range. A linear fit to the data 
(L = 1) yields a slope of the decay rate coefficient of 1.2(3) ■ 10~^ cm'^/s/Kelvin. 

We repeated this measurement for molecules in the lowest hyperfine state | — 4, 3/2) (open 
triangles in Fig.[2]^B)). The data again show /? oc T with a slope of 1.1(3) ■ 10^^ cm'^/s/Kelvin, 
similar to that measured for molecules in the | — 4, 1/2) state. However, in the case of | —4, 3/2) 
molecules, hyperfine state-changing collisions are no longer possible and the only possible loss 
channels are the chemical reactions discussed above. Thus, we find that the rate of chemical 
reactions is determined by the p-wave angular momentum barrier and the chemical reaction 
barrier must be below the collision energy. This suggests that these reactions are barrierless and 
can thus occur freely at ultralow temperatures. Meanwhile, the fact that the same loss rate is 
observed for both | — 4, 1/2) and | — 4, 3/2) state molecules suggests that chemical reactions 
dominate the loss in these ground-state molecular collisions. 

To understand the loss rates, we use two models: one is a simple quantum threshold model 
(QT) and the second is a model that uses the formalism of multi-channel quantum defect theory 
(MQDT). In the QT model, the loss rate for collisions with energy equal to or above the height 
of the j9-wave barrier is determined by the Langevin capture rate (40), which assumes that the 
probability for chemical reactions/hyperfine state-changing collisions is unity. For all collision 
energies below the height of the p-wave barrier, we assume in this model that the loss follows 
the Bethe-Wigner threshold laws (6, 7). Using this assumption, we obtain a simple analytical 
expression for the p-wave loss rate coefficient of two indistinguishable molecules, which scales 
linearly with T (41): 
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where is the reduced mass and h is the Planck's constant divided by 27r. Using a van der 
Waals dispersion coefficient of Cq=\6\30 a. u. for KRb-KRb with a ±10% uncertainty (34, 35), 
the slope of the rate coefficient is predicted to be 1.5(1) • 10~^ cm^/s/Kelvin, which agrees well 
with the experimental measurement. 

In the second model, the loss rate coefficient is found directly by calculating the quan- 
tum tunneling rate through the p-wave barrier (42). This calculation gives (3 — 0.8(1) • 10~^ 
cm^/s/Kelvin, which agrees with the experiment within mutual uncertainties. This P can also 
be derived analytically from the properties of the long-range potential to give (3 = (11.48a)^ 
(ksT/Zi), where a = 0.4778(2/xC6/^^)^/^ = 6.3 nm is the characteristic length of the van der 
Waals potential (11). This fully quantum calculation can be put in the same form as the QT 
model Eq. (2), and gives a (3 that is smaller by a factor of 0.528. With these simple theories, the 
agreement with our molecule-molecule coUisional loss measurements suggests that the chemi- 
cal reaction rates are strongly influenced by the long-range interactions. This observation opens 
intriguing control possibilities since the long-range interaction can be controlled by selecting 
quantum states and tuning collision energies via applied electric and magnetic fields. 

Reaction rates should be dramatically different if molecules are prepared in a mixture of 
different hyperfine states as s-wave scattering becomes allowed. We measured the inelastic col- 
lision rates for ro-vibronic ground-state molecules that were prepared in a roughly 50-50 inco- 
herent mixture of the two hyperfine spin states | — 4, 3/2) and | — 4, 1/2). The time-dependent 
number density of trapped molecules was measured for both spin states. We observed the 
same loss rate for both states, consistent with loss due to collisions between distinguishable 
molecules in different spin states. The rate coefficient is determined to be 1.9(4) • 10~^° cm^/s. 
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independent of temperature, as shown in closed squares in Fig.|2]^B). In comparison to our mea- 
surements for p-wave collisions between spin-polarized molecules, the s-wave collision rate 
between molecules in different hyperfine states is 10 to 100 times larger for a similar tempera- 
ture range. 

The MQDT model can also be used to estimate collision rates where the dominant colli- 
sion channel is s-wave. Here, the relevant length scale is determined by the inter-molecular 
van der Waals potential without any angular momentum barrier (43). We assume that when 
the molecules approach each other within the van der Waals length a, chemical reactions take 
place and remove these entrant molecules with a near-unity probability. The universal loss rate 
coefficient, /3 = 2{h/ ii)a {11), predicts a (3 value of 0.8 ■ 10^^*^ cm'^/s, which is a factor of 2 
lower than the experimentally observed value. This difference suggests that short-range physics 
may play some role in the loss dynamics. An enhancement in the rate coefficient (up to the 
energy-dependent unitarity limit, which is 4 ■ 10^^° cm'^/s for a gas at 400 nK) is possible if 
there is a partial reflection of the colliding species back into the entrance channel (42). The 
reflected amplitude interferes with the incoming amplitude and can either increase or decrease 
the rate coefficient from its "universal" value. Additional theory and experiment are needed to 
explore this possibility. 

In addition to molecule-molecule reactions, the relatively long lifetime of a pure gas of 
spin-polarized molecules in the optical trap (~1 s) affords time to look for chemical reactions 
between atoms and molecules. The MQDT model described above can also predict the atom- 
molecule reaction rates determined by the long-range physics for the universal loss mechanism. 
To prepare the atom-molecule mixture, we control the atom density by selectively removing 
or heating unpaired atoms after the initial molecule creation {44). For these experiments, we 
typically work with an atom number about 5 to 15 times larger than the molecule number. All 
atoms and molecules are prepared in their lowest energy states at 545.9 G. Specifically, K atoms 
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are in their \F = 9/2, nip = -9/2) state, Rb in |F = 1, mj. = 1), and KRb in | - 4, 3/2 >. 
Here F is the total atomic spin and nip is the spin projection. We performed two separate 
experiments, one with K and KRb and the second with Rb and KRb, at a temperature less than 
IfiK. In both cases, the background of atoms in the other, undesired, species was less than 1000, 
corresponding to a density below 5x10^ cm^^. Because both the atoms and the molecules are 
prepared in their lowest energy states, trap loss due to inelastic spin-changing collisions is not 
possible. However, from Table [1} we do expect loss from chemical reactions for the exothermic 
reaction K + KRb K2 + Rb, whereas the endothermic Rb + KRh — > Rb2 + K should be 
forbidden. 

For each experiment, we measured the time dependence of the trapped molecule population. 
Typical molecular loss curves are shown in Fig. |3j^A). To extract the inelastic collision rate, 
we assume that the atom number density is constant. (This is approximately true because the 
number of atoms is much larger than the number of molecules.) The trapped molecule number 
should then decay as 

~77 "^molecule ' -^molecule ' ^atoiru (3) 

at 

where A^moiecuie is the number of molecules, ^atom is the atomic density, and (3 is the inelastic 
rate coefficient. We can then extract /3 via an exponential fit, e~^ "=''°™ *. As can be seen in 
Fig.[3];B), in general we find that the molecules are lost from the trap at a much faster rate when 
K atoms are present than when a similar density of Rb is present. 

To see if the molecule loss arises from atom-molecule collisions, we measure the depen- 
dence of the loss rate on the atom gas density. For the case of KRb -i- K, we observe a clear 
linear dependence on the atomic density and extract an inelastic collision rate coefficient of 
1.7(3) ■ 10~^° cm'^/s. This rate coefficient (due to chemical reactions) can again be predicted 
from the van der Waals length characterizing the long-range part of the potential between the 
collision partners. Using similar methods to those used for KRb -i- KRb, Kotochigova has cal- 
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culated Cq for KRb + K to be 7020(700) a. u. (35), which gives a predicted reaction rate of 
1.1 • 10^^° cm^/s. Again we find that the MQDT value is close but somewhat lower than the 
measured value; this result further solidifies the importance of long-range quantum scattering 
in the inelastic process. 

In contrast, for KRb + Kb collisions, the density dependence of the loss rate is not obvious 
(Fig.[3]^B)). Nevertheless, we again fit the dependence as linear and obtain an upper hmit for the 
rate coefficient of 0.13(4) ■ 10^^° cm^/s, which is one order of magnitude smaller than what we 
measure for KRb + K. 

Our measurement is consistent with the fact that for KRb + Rb, there is no two-body chem- 
ical reaction pathway. A possible mechanism for the residual non-zero rate coefficient for KRb 
-I- Rb would be collisions of ground-state molecules with undetected molecules in high-lying 
vibrationally excited states. These contaminant molecules could be produced by inelastic col- 
lisions of Rb atoms with our weakly bound KRb molecules (45) before the two-photon Raman 
process that produces ground-state molecules. Another possible loss mechanism is Rb -i- Rb -i- 
KRb three -body inelastic collisions. Further experiments will be needed to check these pos- 
sibilities. We note that this suppressed loss rate is observed only when both KRb and Rb are 
prepared in their lowest energy internal states at 545.9 G. If either the molecules or the atoms 
are in an excited hyperfine state, then the observed inelastic rate coefficient rises again to the 
order of 10^^° cm'^/s. Therefore, for practical reasons, removing Rb atoms is important for 
creating a long-lived sample of ground- state KRb molecules. 

Together the studies presented here show that we have observed barrierless chemical reac- 
tions in the short range, with the rates determined by long-range scattering dynamics dictated 
by quantum statistics, angular momentum barriers, and threshold laws. We see that a change as 
seemingly insignificant as flipping a single nuclear spin dramatically changes the rate of molec- 
ular collisions. A pure gas of spin-polarized KRb molecules is long-lived in an optical trap 
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(surviving for a time on the order of a second) whereas a spin mixed KRb sample decays 10 to 
100 times faster. Our results clearly show that chemical reactions can proceed with high rates 
in the ultracold regime. 
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Figure 1: Hyperfine structure of ro-vibronic ground-state ^°K^''Rb molecules at 545.9 G. We 
label the 36 nuclear spin states by their spin projections, mf^ and mf. The energy spacing 
between hyperfine states is ~ h-l30 kHz for |Am^| = 1 and ~ h-760 kHz for \ Amf'''\ = 1. 
By comparison, at a temperature of 300 nK, the molecules' thermal energy is equivalent to 
~ /i-6 kHz, which is more than an order of magnitude smaller than the spin flip energy. In our 
experiments, molecules are prepared in either a single state or in a mixture of | — 4, 1/2) and 
the lowest energy state | — 4, 3/2) (open ellipses). 
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Figure 2: Inelastic collisions between spin-polarized (indistinguishable) or different spin-state 
(distinguishable) fermionic molecules in the ro-vibronic ground state of ^°K'^^Rb. (A) Sample 
data showing the time dependence of the molecule number density. Here the molecules are 
prepared in a single hyperfine state, | — 4, 1/2), and the molecular density decays slowly with a 
rate coefficient of 3.3(7) ■ 10^^^ cm^/s at T = 250 nK. (B) Loss rate coefficient vs temperature. 
The loss rate increases linearly with temperature for spin-polarized molecules, which collide 
via p-wave (inset to A) at low temperature. Data were taken for molecules prepared in either 
I — 4, 1/2) (closed circles) or the lowest energy state | — 4, 3/2) (open triangles). A linear fit 
(solid line) to the | — 4, 1/2) data yields the temperature-dependent loss rate to be 1.2(3)-10~^ 
cm^/s/Kelvin. For the | — 4,3/2) case, where the coUisional loss can only be due to chem- 
ically reactive scattering, the loss rate is similar. The dotted and dashed lines are theoretical 
predictions from the QT model and MQDT, respectively. When the molecules are prepared in 
a mixture of the | — 4, 1/2) and | — 4, 3/2) states (filled squares), we observe a temperature- 
independent decay rate that is 10 to 100 times higher than for the spin-polarized case. The error 
bars represent 1 SD of the decay rate coefficients arising from fluctuations of the molecular 
density, temperature, and fitting uncertainty of the two-body loss curves. 
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Figure 3: Collisions of atoms and molecules in their lowest energy internal states. (A) Sample 
decay curves for the molecular number when subject to collisions with K atoms (circles) and 
with Rb atoms (triangles). The atom numbers are 5 to 15 times higher than the molecule num- 
bers. The reduced initial molecule number for KRb+Rb is due to the coUisional quenching of 
the initial weakly bound KRb molecules with Rb (45) before the two-photon Raman process that 
transfers the molecules down to the ro-vibronic ground state. (B) Dependence of the molecule 
loss rate on the atomic gas density. We observed strong molecule loss with /? = 1.7(3) • 10^^° 
cm^/s for KRb-i-K collisions and suppressed loss of KRb for KRb-i-Rb collisions. The error bars 
are 1 SD of decay rates and atomic densities. 
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